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^  ABSTRACT 

The  polynomial  space  'ijrTn  the  sp>ian  of  the  integer  translates  of  a  box  spline  M  admits  a 
well-known  characterization  as  the  joint  kernel  of  a  set  of  homogeneous  differential  operators  with 
constant  coefficients.  The  dual  space  TC  has  a  convenient  representation  by  a  polynomial  space 
V,  explicitly  known,  which  plays  an  important  role  in  box  spline  theory  as  well  as  in  multivariate 
I)olynomial  interpolation.  '-f  ^ 

7 - iirtliis  paper^we  characterize^the  dual  space  P  as  the  joint  kernel  of  simple  differential  operators, 

each  one  a  power  of  a  direction^  derivative.  Various  applications  of  this  result  to  multivariate 
polynomial  interpolation,  multivariate  splines  and  duality  between  polynomial  and  exponential 
spaces  are  discussed,  ^  ^ 
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On  two  polynomial  spaces  associated  with  a  box  spline 
Carl  de  Boor,  Nira  Dyn  ii.  Amos  Ron 


1.  Introduction 


The  space  E{<f)  of  all  exponentials  in  the  linear  span  S(<^)  of  the  integer  translates  of  a 
compactly  supported  distribution  <j>  is  of  basic  importance  in  multivariate  spline  theory  since,  in 
principle,  it  allows  the  construction  of  good  approximation  maps  to  S{(i>)  from  spaces  containing 
S{<f>).  Generically,  R (0)  is  £)-invariant  (i.e.,  closed  under  differentiation),  hence  is  the  joint  kernel 
for  a  set  {p(D)  :  p  €  of  differential  operators  with  constant  coefficients,  with  ^ 

polynomial  ideal  of  finite  codimension  (in  the  space  x  of  all  multivariate  polynomials,  i.e.,  an 
ideal  of  transcendental  dimension  0,  hence  with  finite  variety).  An  understanding  of  the  interplay 
between  the  space  H(4>)  and  its  associated  ideal  is  useful  in  the  determination  of  the  basic 

properties  of  H{<f))  such  as  its  spectrum,  its  dimension,  and  its  local  approximation  order. 

For  the  important  special  case  when  4>  is  a.  polynomial  box  spline  (and  7i  :=  H{4>)  is  thus  a 
polynomial  space),  an  explicit  set  of  generators  for  the  ideal  is  known  [BHi],  but  nevertheless, 
the  construction  of  their  joint  kernel  was  found  to  be  very  difficult.  At  the  same  time,  a  polynomial 
space  V  (of  very  simple  structure)  is  known  which  serves  as  a  natural  dual  for  7i  and  is  of  substantial 
use  in  the  analysis  of  H.  Specifically,  the  duality  between  H  and  V  has  been  used  in  [DRi]  in  the 
investigation  of  tlie  local  appro.ximation  order  of  some  exponential  spaces,  in  [DR2,2]  in  the  solution 
of  an  interpolation  problem  induced  by  H,  in  [J]  in  the  construction  of  linear  projectors  onto  a  box 
spline  space,  and  in  [DRi]  in  the  computation  of  the  homogeneous  degrees  of  H  (which  is  equivalent 
to  computing  the  Hilbert  function  of  Iji).  See  also  [DM3;  §3]. 

It  is  the  purpose  of  this  paper  to  establish  the  surprising  result  that  V,  too,  is  the  joint  kernel 
of  a  rather  simple  set  of  constant  coefficient  differential  operators,  each  being  just  a  power  of  a 
directional  derivative.  This  allows  us  to  characterize  V  in  terms  of  the  degrees  of  its  polynomials 
when  restricted  to  certain  h'near  manifolds.  Various  applications  of  this  result  to  multivariate 
polynomial  interpolation,  box  spline  theory,  and  duality  between  polynomial  and  exponential  spaces 
are  discussed  as  well. 

In  section  2,  after  defining  the  space  V  and  its  associated  differential  operators,  we  prove  that 
V  is  indeed  the  joint  kernel  of  these  operators.  In  section  3,  we  identify  P  as  a  space  of  least 
degree  among  all  polynomial  spaces  that  interpolate  correctly  on  certain  subsets  of  the  integer 
lattice.  As  a  matter  of  fact,  the  discussion  in  that  section  may  have  an  independent  value:  this 
discussion  illustrates  how  the  interpolating  space  of  least  degree  from  [BRi]  may  be  computed 
using  the  technique  from  [BR2]  of  “perturbing  the  generators  of  a  homogeneous  ideal”,  hence  in  a 


computationally  painless  way. 


Section  4  is  devoted  to  the  more  general  discussion  of  duality  between  finite  dimensional 
polynomial  and  exponential  spaces,  a  discussion  which  improves  proofs  and  results  from  [DM2]  and 
[DRi].  Finally,  we  discuss  in  section  5  the  construction  of  piecewise  polynomials  on  the  support  of 


a  box  spline  and  improve  thereby  an  observation  in  (R). 


-/  Cedes 


(Avail  and/or 
Dial  I  Special 


1 


□  □ 


2.  The  Main  Result 


Let  X  be  a  multiset  of  vectors  in  IR'’\0.  We  will  at  times  think  of  X,  equivalently,  as  a  real 
matrix,  of  size  (s  x  #X).  Let  IH(Jr)  denote  the  collection  of  all  hyperplanes  (i.e.,  linear  subspaces 
of  codimension  1)  which  are  spanned  by  some  columns  of  X.  We  associate  with  each  h  €  IH(A’)  a 
nontrivial  linear  polynomial  which  vanishes  on  h,  and  write  this  polynomial 

thus  using  h-^  to  stand  for  any  particular  nonzero  vector  normal  to  h.  We  are  interested  in  the 
ideal 

generated  by  all  polynomials  of  the  form 

(2.1)  Pk  :=  Pk,x  :=  (/i^-,  h  6  1H(X), 

where  X\h  =  {x  £  X  \  x  ^  h}. 

We  assume  that  X  spans  all  of  IR*.  Then  the  only  point  common  to  all  /i  G  IH(A')  is  0,  and 
consequently,  the  variety  of  (i.e.,  the  set  of  common  zeros  of  all  the  polynomials  in  consists 
of  0  alone: 

V/x  =  {0}. 

This  implies  that  the  codimension  of  in  the  space  tt  of  all  polynomials  in  s  variables  (i.e.,  the 
dimension  of  the  quotient  space  v/I^)  is  finite,  and  that  its  kernel  is  a  finite-dimensional 
polynomial  space,  whose  dimension  equals  the  codimension  of  .  Moreover,  1^1  is  stratified, 
i.e.,  spanned  by  its  homogeneous  elements,  since  is  generated  by  homogeneous  polynomials. 
[BR2]  is  a  ready  reference  for  these  known  facts. 

Here,  to  recall  the  definition,  the  kernel  IL  of  an  ideal  /  is  the  set 

(2.2)  {/  G  P'dR*) :  piD)f  =  0,  Vp  G  /} 

of  all  distributions  annihilated  by  the  set  of  differential  operators  induced  by  7.  In  particular,  since 
is  generated  by 

p,  =  {/i^,.)#(^\'*),  /iGni(x), 

consists  of  the  solutions  /  of  the  system  of  linear  differential  equations 

(2.3)  (Z?hi.  =  0,  Vh  G  IH(J^:). 

This  section  is  devoted  to  a  proof  of  the  fact  that  equals  the  polynomial  space 

(2.4)  V{X)  ;=  span{pv  :  V  C  X,  span(X\V^)  =  HI'’}, 
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with 


pv  := 

«€V' 

In  the  proof,  the  multiset 

(2.5)  :={BCX  :  B  invertible) 

of  all  bases  contained  in  X  plays  an  important  role.  We  use  the  abbreviation 

(2.6)  b{X)-.=  m{X). 

(2.7)  Theorem.  The  kernel  /''j.  of  coincides  with  P(X),  and 

(2.8)  dim  /^1  =  b(X). 


The  theorem  follows  from  the  following  three  lemmata. 


(2.9)  Lemma. 


b(X)<  diinP(X). 


(2.10)  Lemma. 


P(X)  C  /^i. 


(2.11)  Lemma. 


dim  I^±  <  b(X). 


Proof  of  (2.9)Lemma[DRi];  Every  polynomial 

9b:=  n 

iGX\b 

with  B  G  IB(-V)  and  (Aj)  arbitrary  constants,  is  in  P{X),  as  follows  readily  by  multiplying  out. 
For  each  B  G  1B(X),  there  is  a  unique  common  zero  of  the  s  linear  factors  (i,  •)  —  Xx,  x  ^  B,  which 
do  not  occur  in  the  associated  qs',  call  that  point  0b-  Choose  now,  as  we  may,  the  constants  (A^)  in 
such  a  way  that  63  ^  Ob’  whenever  B  B' .  (In  fact,  almost  every  choice  of  the  Ax  would  satisfy 
this  condition.)  It  then  follows  that 

9b(^b')  =  0  B  ^  B', 

proving  the  linear  independence  of  the  collection  (<7s)s6B{X)  P{X)-  ♦ 

Proof  of  (2.10)Lemma:  We  have  to  prove  that,  for  each  h  G  1H(X),  Ph(D)  =  (D/jX 
annihilates  P(-V),  i.e.,  that  pii(D)pv  =  0  for  every  V  C  X  for  which  X\V  contains  a  basis.  For 
this,  we  note  that  Dt^i.pvnh  =  0.  hence 

Ph{D)pv  =  {pvnh)ph{D)pyS^^. 
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On  the  other  hand,  since  X\V  contains  a  basis,  V\h  cannot  coincide  with  X\h,  hence  degp^^v^^  < 
#A'\h  =  degp/i  and  therefore  Ph{D)Py\f^  =0.  ♦ 

Proof  of  (2.11)Lemma:  We  prove  the  lemma  by  induction  on  ifX.  For  the  case  #X  =  s  we  ob¬ 
serve  that  is  generated  by  5  linearly  independent  linear  homogeneous  polynomials,  consequently 
contains  only  constants,  and  so  dim/^±  =  1  =  b{X).  Assume  now  that  ^X  >  s.  We  follow 
the  argument  used  in  the  proof  of  [DM3;  Thm.3.1],  decompose  X  as 

X  =:  X'Vx, 


with  spanX'  =  IR*,  and  consider  the  map 

T  :  — »•  X/,giH(A’')-PA  •  9  iPliiD)Q)h€M(X')i 

■where  p ft  :=  Ph,x‘  =  (h-*-,  are  the  generators  of  and  Ph  :=  pti(D){I^±).  Then 

dim <  dim/^  X -I-  ^  dimP/i, 

AeiHCX') 


since 

kerT  =  n  (7^1)  C  7^'x. 

Consequently,  by  induction, 


This  proves  that 
provided  we  can  prove  that 


dim(/'''l)  <  b{X')  +  E  dim  Ph. 
dim  7^1  <  6(A), 


Y,  dim  Ph<#{Be  1B(A) :  x  e  B}. 

helH(X') 

In  particular,  it  is  sufficient  to  prove  that,  for  all  6  €  IH(A'), 
(2.12)  PhC7^'*X 

with 

Xh  :=  (An/i)Ui. 

For,  (2.12)  implies  that 

dim  Ph  <  dim7'’^'‘±  <  b{Xh), 

(the  last  inequality  by  induction),  while 

6(A,)  =  #{PeIB(A):  i€P}. 

/i€IH(X') 
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The  claim  (2.12)  is  trivial  in  case  x  £  h,  since  then  X'\h  =  X\h,  therefore  pk,x'  =  Ph,x  and 
so  P/,  =  {0}  in  that  case. 

We  now  prove  (2.12)  for  the  contrary  case,  i.e.,  the  case  when  x  ^  h.  We  have  to  show  that 
for  every  €  lH(X;i) 

(2.13)  pk,xAD)PH  =  {0}, 
with 

U  X  ^  k,  then  k  =  h  and  (since  Xh\h  =  {x})  Pk,x^  is  a  linear  polynomial,  and  there  is  nothing 
to  prove  since  then  Pk.x^Ph  X’  =  Ph,Xi  while  pk,x  annihilates  1^1.  by  definition  of  .  For  the 
contrary  case,  i.e.,  the  case  that  x  £  k,  v/e  need  to  prove  that 

(2.14)  P^xP;^x~'"(/^-L)  =  0, 
with 

m:=#(X,\fc),  >l-m:=#(XV). 

For  this,  it  is  sufficient  to  show  that 

(2.15)  g  ideal{I“<^i},, 

with  I  running  over  all  hyperplanes  spanned  by  elements  of  X  and  containing  fc  n  h,  and  with 

K:={k^,.),  P:=(/i^,.),  T:=(f^,-), 

a(0  :=  #(X\f). 

For,  (2.14)  follows  from  (2.15)  since  each  generator  of  the  ideal  in  (2.15)  appears  also  in  the 
defining  set  (2.1)  of  generators  of  ,  hence  annihilates 

We  prove  (2.15)  by  writing  each  polynomial  L  as  linear  combination  of  the  linear  polynomials 
K  and  P,  thereby  obtaining  the  general  homogeneous  element  of  degree  A  of  the  ideal  in  (2.15)  in 
the  form 

Y.iaK  + 

e 

with  Ti  a  homogeneous  polynomial  of  degree 

/i(f)  :=  .4  -  a(f), 

all  i.  We  then  show  that  the  resulting  linear  system  (for  the  coefficients  of  the  various  rt)  has 
a  solution  by  showing  that  its  coefficient  matrix  is  the  transpose  of  the  matrix  which  occurs  in 
(univariate  polynomial)  Hermite  interpolation. 

Here  are  the  details. 
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Since  each  I  contains  k  C\  h,  can  be  written  uniquely  as  a  linear  combination  of  and 
/»-*-.  We  find  it  convenient  in  the  sequel  to  have  the  weights  in  this  linear  combination  sum  to  1, 
i.e.,  to  have  in  the  affine  hull  of  and  h^.  This  we  can  achieve  by  first  choosing  the  (signed) 
magnitudes  of  the  nonzero  vectors  k^  and  so  that  their  difference  is  not  perpendicular  to  any 
of  the  finitely  many  t.  Then,  for  each  £,  we  choose  the  nonzero  vector  to  lie  on  the  line  through 
k^  and  i.e.,  so  that 

=:  {I  -  P)k^  + 

for  some  /?  =  Then  Z  =  (1  —  P)K  +  PH ^  hence 

Z“  = 

j 

with 

B'(t)  :=  ())(1 

the  polynomials  appearing  in  the  Bernstein  form.  Since  DBj  =  we  have  more 

generally 

{(a  +  iy./a\){H  -  Kyi^  =  ^  A'-’^“+‘--’Z)’B“+‘(/3). 

i 

This  means  that  we  have  available  in  our  ideal  an  expression  of  the  form 

M(f) 

ED'f-  A')‘Z^-‘c(f,t)  =  D'Bf{P)c{t,i) 

r  t=o  j  e  i=o 

to  match  the  monomial  H  .  Such  a  match  is  possible  provided  the  linear  system  of  conditions 
imposed  thereby  on  the  coefficients  c(f,i)  is  solvable. 

We  begin  the  proof  that  this  linear  system  is  indeed  solvable  by  showing  that  its  coefficient 
matrix  is  square.  With  the  abbreviation  k'  :=  hO  k,  we  compute 

m  =  #(XAfc)  =  #((^  n  h)\k')  =  #((X'  n  h)\k% 


therefore 

(2.16)  >1  =  #((A"  n  h)\k')  +  #{X'\h)  =  #{X'\k'). 

Also, 

^liC)  =  A-  aiC)  =  #(A"U*')  -  mV)  =  -1  +  m\k')  -  #(A^V)  =  #((^  n  e)\k')  -  1. 

Therefore,  the  number  of  unknowns  is 

+  1)  =  E  ^  =  #(^\^')  =  +  1. 

e  ( 
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i.e.,  equal  to  the  number  of  equations.  Here,  the  sums  are  over  all  £  €  IH(A')  which  contain  k' , 
which  implies  that  X\k'  is  the  disjoint  union  of  the  sets  (X  n  C)\k'  and  so  justifies  the  second 
equality. 

Now  organize  the  unknowns  by  i  and,  within  £,  by  i  =  0, . . .  and  order  the  equations  by 
j  =  0,...,A.  Then  the  matrix  consists,  more  precisely,  of  one  block  of  columns  for  each  £,  with 
the  ith  column  (in  the  block  for  £)  containing  the  value  at  (3  =  Pt  of  the  ith  derivative  of  all  the 
polynomials  Bf,  j  =  0, A,  i  =  0,...,/x(^).  Hence  our  matrix  is  the  transpose  of  the  matrix 
which  occurs  in  the  linear  system  for  the  determination  of  the  Bernstein  form  of  the  polynomial  in 
TT^  which  agrees  with  some  function  (^(^)  +  l)-fold  at  /3  =  Pe,  all  1.  Since  such  univariate  Hermite 
interpolation  is  correct  (since  P^  ^  /3<-  for  I  i'),  the  invertibility  of  our  matrix  follows.  4 

We  note  that  (2.7)Theorem  now  allows  us  to  conclude  that  all  inequabties  appearing  in  its 
proof  must  be  equalities.  This  implies,  e.g.,  that 

jX  ^  jX  whenever  X'  C  X, 

and  that  ph^x'iD)(I^±)  =  Another  immediate  consequence  is  the  following 

(2.17)  Corollary [DM3], [DRi]. 

dimViX)  =  b{X). 

Furthermore, 


(2.18)  d(.V):=  min{#(.Y\/j):  6  ni(A:)} 

is  the  least  degree  of  the  generators  of  7'’^,  hence,  since  =  V{X),  we  have  the  following. 

(2.19)  Corollary [DRi].  With  d{X)  as  in  (2.18), 

T^<d{X)  C  ViX), 


but 

MX)  V{X). 

By  its  definition,  /''±  is  the  set  of  all  polynomials  p  €  tt  that  satisfy  the  foUowing  condition: 

(2.20)  Condition.  For  every  h  €  ni(A')  and  v  €  IR*,  P|v+/ix  C  subspace 

orthogonal  to  h). 

On  the  other  hand,  every  p  €  V{X)  satisfies  an  apparently  stronger  condition  (cf.  [DRoiProp.l]): 

(2.21)  Condition.  For  every  siibspace  M  of  fll*  which  is  spanned  by  elements  from  X ,  and  for 
every  v  €  IR^, 

Plv+/lf±  €  ’’■#(, v\M)-dmi  An- 

Hence  we  conclude  the  following  from  (2.7)Theorem; 
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(2.22)  Corollary.  The  space  P(X)  is  characterized  either  by  (2.20)Condition  or  by  (2.21)Condi- 
tion.  In  particular,  these  two  conditions  are  equivalent  on  tt. 

(2.22)  Corollary  verifies  the  claim  made  in  (DR2;  Remark  2]  that  (2.21)Condition  characterizes 
V(X).  In  ca^e  X  consists  of  N  repetitions  of  5  +  1  vectors  Y  C  iR*  in  general  position,  the 
characterization  of  V{X)  by  (2.21)Condition  has  been  proved  in  [G]  by  other  methods. 

3.  An  associated  polynomial  interpolation  problem 

Here  we  identify  certain  e.xponential  spaces  H  whose  corresponding  “limit  at  the  origin”  H i 
coincides  with  V{X)  (for  an  appropriate  choice  of  X),  and  use  this  identification  in  the  solution  of 
an  associated  interpolation  problem.  The  map 

(3.1) 

which  associates  with  every  finite-dimensional  space  of  entire  functions  a  homogeneous  space  of 
polynomials  of  the  same  dimension,  has  been  introduced  and  studied  in  [BRi]  in  the  conte.xt  of 
a  multivariate  polynomial  interpolation  problem,  and  has  been  discussed  cis  well  in  [BRj]  in  the 
context  of  kernels  of  polynomial  ideals.  To  begin  with,  we  recall  the  definition  of  /f  j  and  review 
some  of  the  results  of  [BRi.i]  needed  here.  Then  we  discuss  a  certain  interpolation  problem  and 
its  relation  to  V{X). 

Given  a  function  f  ^  Q  analytic  at  the  origin,  we  write  its  power  series  expansion  at  the  origin 
in  the  form 

/  =  /o  +  /i  +  /2  +  — , 

where,  for  each  j,  fj  is  a  homogeneous  polynomial  of  degree  j,  and  define  /|  :=  fk  with  k  = 
min{j  :  fj  ^  0};  i.e.,  fi  is  the  first  non-trivial  homogenenous  polynomial  in  the  power  expansion 
of  /.  Using  this  notion  of  the  least  term  /j  of  /,  we  then  define 

(3.2)  Hi  :=  span{/i  :  /  e  H}, 
and  have  [BRj] 

(3.3)  dim  //j  —  dim  H. 

We  say  that  a  pointset  ©  C  IR^  is  correct  for  the  polynomial  space  P  if  the  restriction  map 
P  C®  :  p  P|©  is  invertible.  Equivalently,  ©  is  correct  for  P  if,  for  every  data  {dg)e^Q,  there 
is  exactly  one  p  €  P  for  which  p{9)  =  d$  for  all  0  G  ©.  In  other  words,  interpolation  from  P  at  the 
points  of  ©  is  correct. 

Result  [BRi].  ©  is  correct  for  (expQ)i.  Moreover,  among  all  polynomial  spaces  P  for  which  0 
is  correct,  (exp©)!  is  “of  least  degree”  in  the  sense  that 

(3.4)  dim(P  n  tTj)  <  dim((exp0)i  D  tt^)  ,  Vj. 
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Here, 


(3.5)  expQ  ;=  span{e9}eg0, 

where  eg  :  x  >-*  is  the  exponential  function  (with  frequency  6). 

In  view  of  this  result,  it  is  useful  to  be  able  to  identify  the  “least”  space  (exp0)j  for  given  0. 
This  we  now  do  for  certain  pointsets  0  =  i/x  associated  with  the  box  splines  Mx-  Our  tool  is  the 
following  result 

(3.6)  Result[BR2]. 

=(/hi)J-, 

which  obtains  as  the  kernel  of  the  ideal  generated  by  the  leading  terms  of  the  annihilator  of  H . 
Precisely,  [BR2;(4.3)Theorem  (b)]  provides  the  statement  that,  for  any  polynomial  ideal  I  of  finite 
codimension, 

(/l),  =  (/t)±, 

and  (3.6)Result  is  obtained  by  applying  this  to 

(3.7)  /  =  /h:={p€x:  p{D)H  =  0}, 

which  is  an  ideal  since  //  is  closed  under  differentiation.  In  fact,  with  H  =  exp^,  In  consists  of 
all  polynomials  which  vanish  at  0,  hence  has  finite  codimension  and 

H  =  /Hi., 

(cf.  [BRi;  §3]).  From  /h,  we  obtain  its  homogeneous  counterpart  /nf  as 

(3.8)  /ht  :=  span{pt  ;  p  G /h}, 

where  Pf  is  the  leading  term  of  the  polynomial  p,  namely  the  homogeneous  polynomial  satisfying 

deg(p-  p,)  <  degp. 

The  result  (3.6)  is  of  interest  here  since  it  is  Ccisy  to  identify  elements  of  In  in  ca.se  H  =  exp0: 
If  {pj)  are  linear  homogeneous  polynomials  for  which  the  union  of  the  corresponding  hyperplanes 
{x  €  IR*  :  Pj(x)  =  Cj}  (for  suitable  choices  of  the  constants  Cj)  contains  0,  then  p  ;=  \[j{Pj  -  C;)  6 
/h,  hence  OjPj  ^  we  obtain  enough  of  these  p  to  generate  all  of  then  we  know  by 

(3.6)Result  that  // j  is  the  joint  kernel  of  all  the  corresponding  differential  operators  p{D).  In  fact, 
since  we  know  from  (3.3)  that  dim  //  j  =  dim  //,  we  can  already  reach  this  conclusion  when  we  only 
know  that  the  p  so  constructed  generate  an  ideal  J  of  codimension  <  dim  H . 
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(3.9) Re'  alt  [BRo].  If  the  ideal  J  generated  from  the  leading  terms  of  some  polynomials  in  In 
has  codimension  <  dim  H ,  then  J  =  Ifjy,  therefore 

Hi  =  J±. 

In  our  case,  we  have  identified  (in  (2.7)Theorem)  P(X)  as  the  joint  kernel  of  the  differential 
operators  (with  h  running  over  IH(X)),  hence  are  entitled  to  conclude  that  (exp©);  = 

V{X)  whenever  we  can  find,  for  each  such  h,  constants  Cj^h  so  that  ,  ■)  -  Cj^h)  vanishes 

on  0,  and  know  additionally  that  #0  >  dimP(X). 

Just  such  a  pointset  is,  under  certain  assumptions  on  X ,  provided  by 

(3.10)  ux  :=  {«  €  Z*  :  z  -  ct  =  ^  t^x,  0  <  tj;  <  1,  Vi}, 

xex 

with  z  €  1R*\ U/ieiH(A')(^  +  7Z^].  The  set  -l>x  comprises  the  integer  points  in  the  support  of  the 
shifted  box  spline  Mx{-  +  x)  (cf.  [DMt])- 

Assume  now  that  X  is  unimodular,  i.e.,  the  columns  of  X  are  from  ZZ’\0  and  every  B  G  IB(  A') 
has  determinant  ±1.  For  such  unimodular  X,  the  observations  made  in  [R;§4]  (especially  before 
the  proof  of  Th.  4.1  and  in  the  proof  of  Cor.  4.2)  confirm  the  existence,  for  each  h  G  IH(A’’),  of 
consecutive  integers  Ch,j  so  that  i>x  is  contained  in  the  union  of  the  hyperplanes 

{h^,-)  =  Chj,  y  =  l,...,#A'\/i. 

Indeed,  fixing  h  G  1H(A)  and  taking  the  normal  fi-'-  to  be  a  relatively  prime  integer  vector, 
implies  that  |(/i-‘-,x)|  =  1  for  every  x  G  X\h  since  X  is  unimodular.  By  choosing  the  signs  of  the 
columns  of  A'  appropriately  (which  amounts  to  a  shift  in  A/v)>  we  can  achieve  that  {h-^,x)  >  0  for 
all  X  G  X\h,  hence 

Ch  := 

xe.v\A 

Consequently,  with  z  chosen  so  as  to  satisfy  c*  —  1  <  {h^,z)  <  c^,  I'xix)  must  lie  in  the  union  of 
the  hyperplanes 

=  h  j  ~  0,...,c/,  -  1. 

Moreover,  ^ux  =  dimP(A'),  because  of  (2.17)Corollary  and  the  following. 

(3.11)  Result[DM2].  If  X  is  unimodular,  then 

#<^x  =  b{X). 

This  establishes  the  following  theorem. 

(3.12) Theorem.  If  X  is  unimodular,  then  (e.xp,,^)j  =  P{X).  In  particular,  i>x  is  correct  for 
P(X),  and  PiX)  is  of  least  degree  among  all  polynomial  spaces  for  which  uy  is  correct. 

Note  that  only  the  inequality  #ux  >  b{X)  was  needed  in  the  proof  of  (3.12)Theorem.  As  a 
matter  of  fact,  the  converse  inequality  is  a  consequence  of  the  theorem. 
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4.  The  duality  between  'H{X)  and  V{X) 

The  ideal  and  its  kernel  V{X)  are  intimately  related  to  another  ideal  Ix  and  its  kerne’ 
which  play  a  fundamental  role  in  the  theory  of  box  splines.  In  the  following,  we  review  some 
of  the  basics  about  Ix  and  "HiX)  and  draw  several  connections  between  the  two  settings. 

Let  Ix  be  the  ideal  generated  by  all  polynomials  of  the  form  pv  =  V  C  X 

and  span(Ar\F)  /  IR*.  This  means  that 

(4.1)  /x  :=ideal{p^^^:h€lH(X)}. 

Let  HiX)  denote  the  kernel  of  Ix,  i.e., 

n(X):=Ix±. 

It  is  known  [BHi],  [DMi]  that  H{X)  is  a  finite-dimensional  polynomial  space  and  [DM2]  that 

(4.2)  dim>f(X)  =  6(A'). 


The  spaces  'P(X)  and  'H{X)  are  dual  to  each  other  in  the  following  sense.  Each  p  £  tt  gives 
rise  to  a  linear  functional  p'  on  tt  (and  even  on  a  larger  space  of  smooth  functions),  viz.  the  linear 
functional 

p*  :  q  p{D)q{0). 

This  allows  us  to  consider,  for  any  two  finite-dimensional  linear  polynomial  spaces  Q  and  R,  the 
map 

M  :Q  R’  rp^-p'ifi. 

If  M  is  invertible,  we  say  that  Q  is  dual  to  R  (in  the  sense  that  we  can  then  use  the  elements  of 
Q  in  this  fashion  to  represent  uniquely  the  dual  of  R).  Note  that  the  dual  to  M  carries  R"  =  R 
in  the  same  way  to  Q’,  hence  Q  is  dual  to  /Z  iff  i?  is  dual  to  Q. 

A  necessary  and  sufficient  condition  for  such  duality  is  that  M  be  1-1  and  dim /2  <  dimQ 
(since  then  M  is  necessarily  onto).  In  particular,  if  dimQ  =  dim  JS,  then  such  duality  is  assured 
as  soon  as  we  know  that,  for  every  q  £  Q\Q,  there  is  an  r  G  i?  for  which  g*(r)  ^  0.  By  the  duality 
already  mentioned,  this  is  equivalent  to  having,  for  every  r  £  R\0,  a.  q  £  Q  ioi  which  9*(r)  ^  0. 
For  q  £  Q,  the  linear  functional  q‘  cannot  tell  the  difference  between  /  and  T*/  :=  fo  +  ...  + 
:=  the  power  expansion  of  /  up  to  order  k,  if  k  is  sufficiently  large.  This  allows  us  to  extend 
this  notion  of  duality  to  pairs  Q,R  in  which  fZ  is  a  finite-dimensional  space  of  smooth  functions. 
We  will  eventually  make  use  of  the  following  observation: 

(4.3)  Theorem.  Let  P  be  an  n-dimensional  homogeneous  polynomial  space.  Let  H  he  an  n- 
dimensional  space  of  entire  functions.  If  P  is  dual  to  Hi,  then  P  is  dual  to  H . 

Proof:  For  any  /  £  H\0,  fi  £  H i\0,  hence,  by  assumption,  p*(/i)  7^  0  for  some  p  £  P.  Further, 
since  P  is  spanned  by  homogeneous  polynomials,  we  may  assume  that  this  p  is  homogeneous.  But 
then  p*(/)  =  P*(/i)  5^  0,  showing  that  the  linear  map  II  P"  :  f  {p  ^  p’{D)f{0))  is  1-1, 

hence  P  — ♦  //’  :  p  pj'^  is  onto.  Since  dim  H  =  dim  Hi,  and  dim  Hi  =  dim  P  by  assumption,  the 
theorem  follows.  ♦ 
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We  note  that  the  converse  of  (4.3)Theorem  does  not  hold,  in  general.  For,  it  is  easy  to  make 
up  a  nonhomogeneous  polynomial  space  H  together  with  a  homogeneous  P  dual  to  it,  for  which 
the  conditions  dim(7rj  D  P)  =  dim(3rj  n  ffi),  all  j,  fail  to  hold,  while  these  conditions  are  necessary 
for  P  and  Hi  to  be  dual,  according  to  the  following  proposition  of  use  later. 

(4.4) Proposition.  If  the  homogeneous  po/ynom/a/  spaces  Q  and  R  are  dual  to  each  other,  then 

(4.5)  dim(7ry  n  Q)  =  dim(7rj  n  R) 
for  all  j. 

Proof:  Indeed,  if  (4.5)  is  violated  for  some  (minimal)  j  and,  say,  dim(7rjnQ)  >  dim(7rjnil),  then 
there  exists  a  homogeneous  polynomial  q  €  Q  of  degree  j  for  which  q*  vanishes  on  all  homogeneous 
polynomials  in  R  of  degree  j,  and  hence  vanishes  on  all  of  R,  in  contradiction  to  the  duality  between 
Q  and  R.  4 

With  this,  the  meaning  of  the  following  result  is  clear. 

(4.6)  ResuIt[DM3]*,[DRi].  The  polynomial  spaces  V{X)  and  'H(X)  are  dual  to  each  other. 

In  (3.12)Theorem,  the  space  V(X)  has  been  identified  as  the  least  space  for  certain  inter¬ 
polation  problems.  In  [BR2]  the  space  7i(X)  has  been  identified  as  the  least  space  for  other 
interpolation  problems.  We  now  make  use  of  the  duality  between  ViX)  and  HiX)  to  connect 
'H{X)  with  the  interpolation  problems  associated  with  V{X)  and  vice  versa.  As  a  preparation,  we 
procure  a  class  of  spaces  whose  corresponding  least  space  is  'H{X)  in  much  the  same  way  in  which 
we  obtained  suitable  exponential  spaces  exp^^  whose  leaist  is  V{X):  We  perturb  the  linear  factors 
of  the  set  of  generators  for  the  ideal  (viz.  ly)  whose  kernel  is  H{X).  Specifically,  given  the  set 
r  =  {(i,A;r)  \  X  £  X  ,  Ax  €  C},  we  consider  the  ideal 

(4.7)  Ir  :=ideal{qh:h€m{X)}, 
with 

q,:=  n  ((x,.)-Ax). 

x^X\h 

Then  W(r)  :=  /fi.  is  an  exponential  space  (i.e.,  a  space  which  is  spanned  by  certain  products  of 
exponentials  with  polynomials),  and  [BRi] 

(4.8)  W(r),  =  n{X). 

The  diagram  below  illustrates  for  a  unimodular  X  the  various  connections  established  so  far 
between  the  ideals  Ix,  their  kernels  and  the  associated  exponential  spaces. 

*  The  authors  in  [DM3]  attribute  the  result  to  Hakopian. 
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(exp.Ji  =  V{X)  = 

T 

dual 

i 

H(r)i  =  n{x)  =  ixL 

Diagram 


Our  first  corollary  improves  [DM2;Thm.  4.1]; 

(4.9)  Corollary.  Let  X  be  unimodular.  Then  t/x  is  correct  for7i(X),  and  'H{X)  is  of  least  degree 
among  all  polynomials  spaces  for  which  t/x  is  correct. 

Proof:  We  apply  (4.3)Theorem  with  H  =  exp^^  and  P  =  'H{X).  By  (3.12)Theorem  =  V{X), 
while  by  (4.6)Result,  P(X)  is  dual  to  'H(X).  Since  also  H{X)  is  homogeneous  (as  the  kernel  of 
a  homogeneous  ideal),  (4.3)Theorem  implies  that  TL(X)  is  dual  to  exp^^,  which  is  equivalent 
(cf.  [BRi;§4])  to  the  correctness  of  ux  for  7{{X).  Moreover,  since  'P{X)  and  'H{X)  are  both 
homogeneous  and,  by  (4.6)Result,  dual  to  each  other,  we  must  have 

(4.10)  dim(ffj  n  'H{X))  =  dim(xj  n  V{X)) 

for  every  j,  by  (4.4)Proposition.  Since  we  already  know  by  (3.12)Theorem  that  V{X)  is  of  least 
degree  among  all  polynomial  spaces  for  which  i>x  is  correct,  it  follows  from  (4.10)  that  H{X)  has 
the  same  property.  4|k 

(4.11)  Corollary [DRi].  The  spaces  V(X)  and  7i(r)  are  dual  to  each  other. 

Proof:  Take  F  =  P(X)  and  H  =  K(r)  in  (4.3)Theorem.  Since,  by  (4.8),  'H(r)|  =  'H(A'),  and, 
by  (4.6)Result,  7{{X)  is  dual  to  V{X),  (4.3)Theorem  provides  the  desired  result.  4^ 

We  refer  to  [DRi;  §7]  for  a  discussion  of  the  interpolation  conditions  induced  by  7f(r). 

5.  Application  to  box  splines 

In  this  section  we  point  out  some  connections  between  the  results  of  the  previous  sections  and 
the  theory  of  multivariate  splines.  In  the  discussion  here  the  (polynomial)  box  spline  Mx  associated 
with  a  set  of  directions  X  plays  a  central  role.  For  our  purposes,  it  is  sufficient  to  note  that  Mx  is 
a  piecewise-polynomial  function  supported  on 

Qx  :=  {Xt:  <€[0,1]#^} 

and  satisfies 

7r(Mx)  =  H(A’), 

where,  for  a  general  compactly  supported  4>,  the  notation  ^(d*)  stands  for  the  space  of  polynomials 
spanned  by  the  integer  translates  of  <f). 

We  make  use  of  the  following  result,  which  is  a  special  case  of  [R;  Thm.l.lj: 
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(5.1)  Result.  Let  F  be  a  translation-invariant  space  of  polynomials,  and  fi  a  compact  subset  of 
IR*  with  boundary  dQ.  Then  the  following  conditions  are  equivalent: 

(a)  There  exists  a  function  <p  supported  in  fi  and  satisfying  Tr{<())  =  P,  0. 

(b)  For  every  z  G  ni*\  Uae2z*  oi  +  dfl,  the  set 

'■=  t'fi(-j)  :=  {o:  €  ZZ*  :  z  —  a  £  Q} 


is  total  for  P,  i.e.,  no  element  of  P\0  vanishes  on  this  set. 

Note  that  when  we  take  fi  to  be  ilx  =  suppMx,  the  sets  are  identical  with  the  sets  vx 
from  (3.10).  Thus,  by  appealing  to  (3.12)Theorem,  we  deduce  from  the  impUcation  (b)  =>  (a)  of 

(5.1) Result  the  following 

(5.2)  Corollary.  Let  X  be  a  unimodular  set  of  vectors.  Then  there  exists  a  function  (actually 
many)  ijj  =  V’x  which  is  supported  in  Qx  satisfies 

;r(V>)  =  P(X). 

The  above  corollary  provides  no  information  about  the  smoothness  of  the  compactly  supported 
V’.  Yet,  it  is  known  [BH2],  that  at  least  for  the  special  case  of  the  three-direction  mesh  with  equal 
multiplicities,  no  piecewise-P(Ar)  function  supported  on  can  match  the  smoothness  of  the 
corresponding  box  spline  Mx- 

The  box  spline  Mx  is  a  smooth  function  supported  on  Qx-  Hence,  one  may  hope  that  there 
exist  functions  supported  on  Qx  which  are  less  smooth  than  Mx,  yet  their  corresponding  x{<p)  is 
“better”  in  the  sense  that  it  contains  some  of  the  polynomials  of  lower  degrees  which  were  missing 
in  x{Mx)  =  'H{X).  [R;  Cor.4.2]  gives  a  partial  negative  answer  to  that  hope  by  showing  that  for 
a  unimodular  X  and  a  function  <t>  supported  in  Qx,  if  4>{0)  5^  0  and  xj  C  x(<p)  for  some  j,  then 
Xj  C  'H{X)  =  x{Mx)-  The  following  result  improves  that  corollary. 

(5.3)  Corollary.  Let  X  be  an  unimodular  set  of  directions,  Mx  the  corresponding  box  spline. 
Let  <f>  be  a  compactly  supported  function  satisfying 

supp  4>  C  supp  Mx , 


and 

m  0. 


Then,  for  each  j. 


dim(7rj  D  x{d>))  <  dim(7ry  n  x{Mx))- 


Proof:  Let  u  be  one  of  the  sets  t'x(^)  associated  with  X.  By  (3.12)Theorem,  (exp,,)^  =  P(X). 
Now  we  may  apply  (5.1)Result  to  conclude  that  v  is  total  for  x{(p),  and  hence  x(<f>)  can  be  extended 
to  a  space  Q  for  which  1/  is  correct.  Since  'P(X)  is  the  least  space  of  exp,,,  it  satisfies  the  least 
degree  property  (3.4),  thus  we  conclude  that  for  every  j 

(5.4)  dim(7rj  D  x{(i>))  <  dim(jrj  D  Q)  <  dim(xj  fl  P(X)). 
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On  the  other  hand,  by  (4.4)  and  (4.6), 


dim(7rj  n  P(X))  =  dim{7r_,-  D  W(X)),  Vj, 

and  (5.3)  now  follows  from  the  fact  that  ir(Mx)  =  'H(X).  4 

We  note  that  the  result  is  no  longer  valid  if  we  drop  the  unimodularity  assumption  (cf.  [R;  Ex. 
4.1]). 


6.  A  remarkon  (2.11)Lemma 

A  careful  examination  of  the  proof  of  (2.11)Lemma  shows  that  the  details  of  the  connection 
between  the  ideal  and  its  kernel  I^A.  enter  into  the  argument  in  only  a  minor  way.  The  only 
facts  used  are:  (i)  the  map  rr  — *■  L{t)  :  p  p(D)  is  a  ring  homomorphism;  and  (ii)  for  any  basis  B 
in  X,  dim/®±  =  1. 

This  suggests  the  following  result. 

(6.1)  Proposition.  Let  M  :  k  -*  L(V)  be  a  ring-homomorphism  into  the  ring  of  linear  maps  on 
the  linear  space  V.  For  a  given  multiset  X  of  directions,  define 

••=  Hpg/x  kerM(p). 


Then 


^  dim/f/1. 
BelB(X) 


The  proof  of  the  proposition  follows  entirely  that  of  (2.11)Lemma,  with  the  obvious  modifica¬ 
tions  whenever  the  induction  hypothesis  is  applied. 

It  would  be  nice  to  identify  other  settings  rather  than  the  one  utilized  in  this  paper,  where  the 
above  proposition  is  of  use. 
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